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The Source Galerkin Method is a new numerical technique that is being developed to solve Quantum Field 
Theories on the continuum. It is not based on Monte Carlo techniques and has a measure to evaluate relative 
errors. It promises to increase the accuracy and speed of calculations, and takes full advantage of symmetries of 
the theory. The application of this method to the non-linear a model is outlined. 



1. MOTIVATION 

Conventionally, in numerical approaches to 
field theories, integrals for generating functionals 
are evaluated on lattices using Monte Carlo meth- 
ods. These are quite successful but there are still 
constraints due to lack of computer power and 
algorithmic limitations. These limitations be- 
come important for theories with significant con- 
tributions from fermionic interactions beyond the 
quenched approximation, and for actions that are 
not positive definite. Source Galerkin is an alter- 
native computational method that avoids these 
particular problems ||l|,||. It works both on the 
lattice and the continuum and handles fermions 
as easily as bosons Computations require 

significantly less time than Monte Carlo based 
methods. The basic idea of the method is to find 
solutions to the functional differential equations 
for a field theory in the presence of sources. The 
solutions are generating functionals of the theory. 



expansion over the external sources of the theory. 



Z* = exp [ j jxGxyjy + J ji^jxH^xy 



zjyjz + ' ' ' . 



Z* does not satisfy the source differential equa- 
tions of the field theory exactly, but we can re- 
quire that specific functional projections of these 
equations are satisfied. These projections are 
formed by multiplying by test functions in sources 
and integrating this product in source space us- 
ing a fixed measure. Introducing the appropri- 
ate number of test functions, we can solve for 
the parameters of the resulting expressions. This 
method causes the error in assumed solutions to 
approach zero in the mean and so the ansatz sat- 
isfies the equations approximately. The details of 
this will be made clear in what follows. This ap- 
proach belongs to the general class of techniques 
called Galerkin methods 



1.1. Overview of the Method 

To analyze a theory, we start with a Lagrangian 
in the presence of external sources. The source 
functional differential equations satisfied by Z are 
generated in the usual way. The Source Galerkin 
technique consists of a systematic iterative way 
to solve these equations to increasing accuracy. 
We proceed by introducing an ansatz which is an 



2. APPLICATIONS 

The method has been applied to various field 
theories on the lattice and, more recently, on the 
continuum A(/)^, Gross- Neveu, and other two di- 
mensional fermionic models have been studied. 
Here, we discuss an application of the method to 
the 0(3) Nonlinear a Model with d = 2. Since 
this model is asymptotically free, it serves as a 
toy model for our initial approach to Non-Abelian 
gauge theories. 
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2.1. 0(3) Nonlinear a Model {d = 2) 

This model consists of three real fields coupled 
to each other. The Lagrangian is: 



9' 



x{x)S{x) 



(j) and X si'fi the canonical and Lagrange multi- 
plier fields respectively, j and S are the external 
sources for these fields while g is the coupling con- 
stant of the theory. The Schwinger-Dyson equa- 
tions are 



d 



25Z[],S] 5^Z[j,S] 



Sj^ix) 6j'^{x)SS{x) 



-rix)Z[j,S]^0, 



mx)S.nx) y.^iJ^S]-2Six)Z[j,S]^0. 

These are the functional differential equations to 
which we apply the Galerkin procedure. The sim- 
plest ansatz is 



exp 



r{x)G'^\x,y)j''iy) + xo / Six) 



The inner product that will be used for projecting 
residuals on the test functions is 

The subscript denotes the j space, e weights the 
inner product such that higher order terms con- 
tribute less. For this model, the projections are 

(l,(j''(a;'),i?a),)s-0, 

(l,(l,i?e),)s=0. 

Ra and Rc are the residuals for the equations of 
motion and constraint respectively. These equa- 
tions, are solved for Green's functions and La- 
grange multiplier, xoi and are defined in source 
and coordinate space. The space-time depen- 
dence of the Green's functions is approximated 
by truncating a Sine expansion, 

N 

G{x) = G{kh)S{k,h){x). 

k=-N 



where the Sine function is defined as 

Sine functions have many algebraic properties 
that make them very useful to work with 
While this is a very direct approach, this inter- 
polative approximation requires that every extra 
dimension introduces another sum to the prob- 
lem. Using this approximation 0{N'^) num- 
bers are required to store a two point function. 
The memory limitations become a concern when 
OiN^"^) numbers are to be stored for the Jacobian 
with four-point functions. This leads us to seek 
an alternative method which uses general infor- 
mation about the spectral structure of field the- 
ories. Any exact two point function can be rep- 
resented as a sum over free two point functions 
using spectral representations. Consequently, we 
can shift our basis for numerical solutions to reg- 
ulated free Euclidean Green's functions: 

'^(^) TTTd / — 2~ 2~ 

A term in the spectral representations has the 
Sine expansion: 



G(x) « Go 



N 

E 



c{k) exp 



4(6^='' + 1/A2) 



c{k) 



1 



^kh 



1/A2 



d/2 



exp 



The Galerkin equations are solved approximately 
using only one term of the spectral expansion 
with the parameters m* and Gq. The (3 func- 
tion (Figure. |l|) is calculated and compared to 
results. The ansatz is improved. 



Z* = Z* exp 



S{x)D(x~y)S{y) 



+ \ I r{x)j\y)F'^\x-z,y-z)S{z) 



Zs is the simplest ansatz. and D are the ver- 
tex term and the two-point function for the x field 
respectively. Results get close to exact value dra- 
matically with improved ansatz. We have shown 
this in the model by introducing a four-point 
fimction (Figure. 0). 
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3. CONCLUSIONS 




Figure 1. P 
1/N result. 



function vs. coupling, compared to 




Coupling, A 

Figure 2. Comparison of improved ansatz with 
lowest order and exact answer for X(j)'^, d = 1 



The Source Galerkin method promises to ex- 
tend the range of field theoretic problems that 
can be solved with current computers. The low 
order results of the method applied to the non- 
linear a model show asymptotic freedom and are 
close to the large-iV results. The contributions 
from higher order corrections are being studied. 
Understanding this model is of primary impor- 
tance, as it is a step to supersymmetric and 
matrix models. One of the major advantages 
of our method lies in its approach to fermionic 
fields. Fermions are treated identically to bosons 
except for the anti-commutivity of the associ- 
ated sources. Consequently, expansions involving 
fermions and bosons are symmetric in both types 
of fields and there is no analogue to the fermion 
determinant problem. The preliminary results of 
fermionic theories are very promising ||^. Source 
Galerkin takes advantage of the symmetries of a 
theory. The calculations are done in a continuum 
formulation, thus avoiding lattice-specific prob- 
lems. While studying various ways of approx- 
imating Green's functions an interesting appli- 
cation has been investigated. Using Sinc-based 
spectral representations, we can calculate Feyn- 
man diagrams in a perturbation theory. Com- 
plicated diagrams can be evaluated much faster 
than with Monte-Carlo methods, and to higher 
accuracies |8j. 
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